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Abstract 

We derive analytical expressions for the excitation energy of the isoscalar giant 
monopole and quadrupole resonances in finite nuclei, by using the scaling method and 
the extended Thomas-Fermi approach to relativistic mean field theory. We study the 
ability of several non-linear a — to parameter sets of common use in reproducing the 
experimental data. For monopole oscillations the calculations agree better with exper- 
iment when the nuclear matter incompressibility of the relativistic interaction lies in 
the range 220-260 MeV. The breathing-mode energies of the scaling method compare 
satisfactorily with those obtained in relativistic RPA and time-dependent mean field 
calculations. For quadrupole oscillations all the analyzed non-linear parameter sets 
reproduce the empirical trends reasonably well. 
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1 Introduction 



The relativistic mean field (RMF) approach to Quantum Hadrodynamics |l| has become 
a very useful tool for describing ground-state properties of nuclei along the periodic table. 
The simplest model, the linear a — uj model of Walecka ||, describes the nuclear force 
in terms of the exchange of a and uj mesons. It is known that the value of the nuclear 
matter incompressibility is unreasonably high in this linear model (K^ ~ 550 MeV), which 
is a serious drawback for a precise description of some properties of finite nuclei and of 
collective excitations such as the breathing mode (isoscalar giant monopole resonance). The 
problem can be cured by introducing cubic and quartic self-interactions of the a meson || , 
which in particular have the effect of lowering the incompressibility, and the model can be 
refined by adding an isovector p meson. Current non-linear parameter sets, such as the NL3 
parametrization give ground-state binding energies and densities in very good agreement 
with the experimental data, not only for magic nuclei but also for deformed nuclei as well 
as for nuclei far from the stability line. 

The RMF model has also been applied to describe dynamical collective motions in nuclei. 
The basic theory of vibrational states in nuclei, the random-phase approximation (RPA) 
m , has been generalized to the relativistic domain (RRPA) and it has been used in 

calculations of isoscalar giant resonances, to obtain response functions and mean energies 
for several magic nuclei. Small-amplitude collective motions such as the isovector dipole 



oscillation and the isoscalar and isovector quadrupole oscillations [ 10 1 , as well as the isoscalar 



and isovector monopole oscillations [II], have been studied in the time-dependent RMF 
approach. Another approach is based on constrained RMF calculations. It has been applied 



to obtain breathing-mode energies and incompressibilities in the linear [12,13] and non-linear 



14,15) a — uj models. The generator coordinate method, with generating functions that 



are solutions of constrained RMF calculations, has been employed to compute excitation 



energies and transition densities of giant monopole states []TT],p~4| . Other calculations of 
breathing-mode energies in the relativistic framework, see Refs. [^,|I^JT^], have relied on 
the scaling model in combination with the leptodermous expansion of the finite nucleus 



incompressibility derived by Blaizot |18 
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In the non-relativistic framework it is well established that the RPA is the small ampli- 
tude limit of the time-dependent Hartree-Fock approach ||[I9|. In the relativistic case the 
RPA configuration space must include negative energy states from the Dirac sea in order 
to reproduce the results of time-dependent RMF or constrained RMF calculations PJ2D|. 
Paraphrasing the statement, the RRPA corresponds to the small amplitude limit of the 
time-dependent RMF theory in the no-sea approximation when the RRPA includes both 
positive energy particle-hole pairs, and pairs formed from the empty Dirac sea states and 
the occupied Fermi sea states. 

Semiclassical methods in nuclear physics, like the Thomas-Fermi theory, have proven to 
be very helpful for dealing with nuclear properties of global character that vary smoothly 
with the particle number A (e.g., binding energies, densities and their moments) |6| j2T|j2^ . |23| . 
The success of these methods stems from the fact that the shell corrections (quantal effects) 
are small as compared to the smooth part given by the semiclassical calculation. Semiclas- 
sical techniques like nuclear fluid dynamics [^4| and the extended Thomas-Fermi method 
[p5p6| , prr| have been applied to study giant resonances in non-relativistic models. In the 
relativistic context, the nuclear fluid dynamics approach has been utilized, e.g., in Refs. 



28, 29, 30 1 . The authors of Refs. |31|,|32| resorted to a local Lorentz boost and the scaling 



method to study isoscalar giant monopole and quadrupole states in the linear a — uj model. 
The investigations were carried out for nuclear matter (where a Thomas-Fermi approxi- 
mation is exact) |H[ and for symmetric, uncharged finite nuclei [3^] whose densities were 
solved in the relativistic Thomas-Fermi (RTF) approximation. 



The relativistic extended Thomas-Fermi (RETF) method [ p^j3^j35[| is a refinement of 
the RTF method, which incorporates gradient corrections of order h 2 to the pure RTF ap- 
proximation. It was derived only a few years ago and it has since been applied in calculations 
of ground-state binding energies and radii of finite nuclei |]3"l| , |3"t3| , |3"T| and in investigations of 
nuclear surface properties [ ^ , P^J38|J39[| . In the present work we shall use the RTF and 
RETF approaches to calculate the excitation energies of the isoscalar giant monopole and 
quadrupole resonances in spherical nuclei. This will be done by means of the scaling method, 
within the framework of the non-linear a — uj model and the RMF theory. We shall also 
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perform constrained calculations for the monopole state. 

Recently, the basic theory derived in the RTF approach has been applied to discuss the 
virial theorem and to study the breathing-mode energy within the RMF theory [iD|. In 
the present contribution we analyze our self-consistent method in depth. To our knowledge, 
for realistic non-linear parameter sets of the RMF theory, these are the first calculations 
of isoscalar giant resonances in finite nuclei carried out with the scaling method which are 
fully self-consistent (i.e., we do not make use of a leptodermous expansion as in previous 



scaling approaches ||ll^]lq , |T7| ) . Owing to the meson-exchange nature of the relativistic model 
one has to deal with finite range forces, which renders the scaling method more involved 
than, e.g., for non-relativistic zero-range Skyrme forces. Moreover, in contrast to the non- 
relativistic situation, there exist two different densities, namely the baryon and the scalar 
density, in accordance with the fact that one has two types of fields, the vector and the 
scalar field. 

The article is organized as follows. After the introductory remarks, we collect the basic 
expressions of the energy density and the variational equations of the RTF and RETF models 
in Section 2. The third and fourth sections are devoted to the derivation of the equations 
and the discussion of the numerical applications for the giant monopole and quadrupole 
resonances, respectively. The conclusions are laid in the last section. Some technicalities 
and a derivation of the virial theorem for the relativistic model are given in the appendices. 
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2 Energy density and variational equations 



The mean field Hartree energy density of a finite nucleus in the non-linear o — uj model 
reads §g§ 



H = Y t ip\ -ia • V + (3m* + g v V + ^g p Rr 3 + ^eA(l + r 3 ) 



<p% + H{. 



(2.1) 



The relativistic effective mass (or Dirac mass) is defined by m* = m — g s (j), r 3 is the third 
component of the isospin operator, and the subindex i runs over the occupied states (fi of 
the positive energy spectrum. TC{ stands for the free contribution of the meson fields <p, V 
and R associated with the a, u and p mesons, respectively, and of the Coulomb field A: 



H f 



;V0) 2 + m s 2 2 



3 r 4 r 2 



; W)" + m 2 v V 2 



(VR? + mlR 2 



(2.2) 



It is understood that the densities and fields are local quantities that depend on position, 
even if we do not make it explicit in most of our expressions. Units are h — c—1. 

The semiclassical representation of the energy density ( |2.1| ) has a similar structure, 
except that the nucleon variables are the neutron and proton densities (p n and p p ) instead 
of the wave functions. In the RETF approach it reads p4] , |36| , p7| , |38| , p9| 



H = £ + g v Vp + g p Rp 3 + eAp p + Hf 



(2.3) 



where p = p p + p n is the baryon density, p 3 = \{p v — p n ) is the isovector density, and the 
nucleon energy density £ is written as £ = £ + £ 2 with 



£ = E' 



>7T 



&Fg4g + k Fq e Fq ~ ™* 4 In 



*4 i ^Fg + e Fq 



(2.4) 



and 



£ 2 = E [Bi q (k Fq ,m*)(Vp q ) 2 + B 2q {k Fq ,m*) (Vp q ■ Vm*) + B 3q (k Fq , m*)(Vm*) 2 ] . (2.5) 
g 

For each kind of nucleon (q = n, p), the local Fermi momentum k Fq and e Fq are defined by 



k Fq = (37r p q ) 



2 n W3 



eFq = Jk Fq + m* 2 . 



(2.6) 



The coefficients B iq are the following functions of k? q and m* p4| , |39| : 



B 



7T 



lq 



vrr 



B 



m 



2q 



6kp q e 2 



In 



Fg 



rrr 



B- 



3q 



k 2 

^Fq 



247r 2 e| 9 



fc Fq V k FqJ m * 



(2.7) 



The RTF approximation is obtained by neglecting £2 in Eq. ( |2.3| ). The gradients contained 
in £ 2 arise from the RETF corrections of order h 2 to the functional Eq. Naturally, these 
corrections are more important in the nuclear surface region where the densities and the 
fields change more rapidly 

The semiclassical ground-state densities and meson fields are obtained by solving the 
Euler-Lagrange equations 5H/5p q = \i q (with \i q being the chemical potential) coupled to 
the field equations 



(A - m 2 s )<f) 
(A-m 2 v )V 
(A - m 2 p )R 
AA 



-g s p s + H 2 + c0 3 , 



-9vp, 

-g P p3, 

-epp. 



(2.8) 
(2.9) 
(2.10) 
(2.11) 



The semiclassical scalar density in ( |2.8j ) is given by 

5£ , 5£ 2 

+ = PsO + Ps2 



5m* 5m 
m 



E 



2tt 2 



, *2 , faq + ^Fq 



E 



9 L 



B2 ,Ap, + 2 B3 ,A m . + (^-§f)(V,,) 



3g 



+ 2 _ii ( v P ,.v™«) + 1 ^ 



Vm 



*\2 



(2.12) 
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Parenthetically, we would like to mention that the densities p, p s and 8 above are the 
semiclassical counterparts of the quantal densities p = YliYiVh Ps = J^i^PiPVi an d £ = 
E^lh'a-V + f3m*]<p t . 

Since the energy density H is to be integrated over the space to compute the total 
energy, the field equations (|2.8| )- (|2.11| ) can be used to rewrite TC{, e.g., by transforming 
ICVV) 2 + m 2 7 V 2 } into V(— A + m^V = g v Vp (valid, of course, under an integral sign). This 
way, on defining an effective scalar density by 

g s pf = 9sPs - b4> 2 - c0 3 , (2.13) 

TC can be recast as 

U = £ + l -g4pf + \b<P* + \c<P 4 + \gyV P + \gpRp3 + \eA Pv . (2.14) 

This form of TC will be more convenient for facilitating the calculations to be presented 
below. 
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3 Giant monopole resonance 



As far as the giant resonances are dynamical processes one must first describe the nucleus 
from a moving frame [[31 ,[32]. This is a rather technical matter for our present purposes and 
it is left for Appendix A, where we derive the expression of the energy of a nucleus within 
the relativistic model in a frame moving with velocity —v. After performing the scaling 
of the energy in this frame, one obtains general expressions for the two main ingredients 
required for the calculation of the excitation energy of the giant resonance, namely, the 
restoring force and the mass or inertia parameter (Eqs. ( [A. 12 ) and ( A.13|) , respectively). 

The present section proceeds as follows. We begin by introducing our scaling approach 
for the monopole vibration. Next we obtain analytical expressions for the restoring force and 
the mass parameter of the monopole state. The calculational details of the derivatives of the 
meson fields with respect to the collective coordinate of the monopole vibration are reserved 
for Appendix B, where we also discuss the virial theorem (stationarity of the scaled energy) 
for the relativistic model. Next in the section we give a brief summary of the constrained 
approach to the breathing mode, for comparison with the scaling approach. The section 
closes with the discussion of the results of our numerical calculations. 



3.1 Scaling 

Denoting by A the collective coordinate associated with the monopole vibration, a normal- 
ized scaled version of the baryon density is 

p x (r) = A 3 p(Ar). (3.1) 

Accordingly, the local Fermi momentum changes as 

k FqX (r) = [37r%A(r)] 1/3 = \k Fq (\r). (3.2) 

The meson fields <f), V and R and the Coulomb field A are also modified by the scaling due 
to the self-consistent equations ( pT8|) ( pTTT| ) , which will relate the scaled fields to the scaled 
densities. Unfortunately, the meson fields do not scale according to simple laws like ( |3.1| ) 
and (|3.2j ) because of the finite range of the meson interactions. This is most apparent for 
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the scalar field since the scalar density in the source term of Eq. ( |2.8| ) transforms not only 
due to the scaling of k Fq but also of (j) itself (or m*), see Eq. ( |2.12| ) for p B . For reasons that 
will become clear immediately, we shall write the scaled effective mass m x (r) = m — g s (j)\(r) 
in the form 

m* x (r) = \rh*{\r). (3.3) 

Note that the quantity m* carries an implicit dependence on A apart from the parametric 
dependence on Xr. 

On account of Eqs. (|3.2|) and ( j3.3|) the scaled semiclassical energy density E\ = £o\ + £ 2 x 
and scalar density p s \ = p s oA + p S 2X read 

£ x {r) = \%[k Fq {\r),m*(\r)\ + X%[k Fq (Xr), m*(Ar)] = X 4 £(Xr), (3.4) 
Psx(r) = X 3 p s0 [k Fq (Xr),m*(Xr)} + X 3 p s2 [k Fq (Xr),m*(Xr)] = X 3 p s (\r). (3.5) 

The tilded quantities S and p s are given by Eqs. ( |2.4|) , (|2.5| ) and ( p.!2j ) after replacing m* 
by m*. Note the usefulness of ( |3.3| ) to be able to write ( |3.4| j and ( |3.5| ) in this compact form. 
For the scaled total energy density TC\ we obtain 



n x = x 

with the definition 



X£ + ^g s (j)xpf + + + \gvVxp + ^g P Rxp3 + \eA x p P 



(3.6) 



g s pf = g s p s - - ±<fc. (3.7) 



Observe that the same expression is valid regardless of performing the calculations in the 
RETF model or in the RTF model, as the corrections of order h 2 (£2 and p s2 ) scale in the 
same manner as the the Thomas-Fermi terms (£q and p s o). 

3.2 Restoring force 

The restoring force Cm of the monopole vibration is given by the second derivative of the 
scaled energy with respect to the collective coordinate A, calculated at A = 1 (Appendix A). 
The first derivative of the scaled energy is 



= 1 ^a<i>\ 1 . dpf b 3 



3 c 
4A 1 ' 



a + 



dV x 1 dR x 1 dA x 



(3. 



Here we have used dS/dX = p s dm*/d\ (as p s = 5£/5m*) and 



dmX 



m* + X- 



dfh* 



-ft 



oq>\ 



(3.9) 



from the definition 
A = 1 we have 



dX "'" ' " dX 
of m* . Differentiating again with respect to A and then setting 



M 



dX 2 



d(Xr 



A 3 



J A=l 



/ 



dr 

3 

A 1 



~ m 9A 



1 ~ eS d 2 <Px 



2 9sPs dX 



■ 2 + ^g s <P\ 



B 2 n ett b r 



d 2 v> 



d 2 A 



dX 



2 9vp dX 2 



a ^ 1 &R X , 1 



(3.10) 



A=l 



The calculation of the derivatives of the scaled meson fields with respect to A is illustrated 
in Appendix B. There we also work out Eq. ( ^.8[ ) at A = 1, which leads to the virial theorem 
(stationarity condition of the energy) for the relativistic mean field model. Following the 
technique outlined in Appendix B we find 



d 2 Vx(r) 



dX 2 



A=l 



dr'p(r') 



2s 



<WM , 2 d 2 V UJ (s) 



ds 



+ s 



ds 2 



with 



VJs) 



9v e" 



(3.11) 



(3.12) 



47T S 

An equivalent expression holds for the field R\. As is well known, the second derivative 
of the scaled Coulomb field A\ vanishes [HJ (you only have to evaluate ( 3.1 1| ) for a zero 
meson mass). For the scalar field one gets a lengthier expression owing to the extra implicit 
dependence of p~l s on A: 



dX 2 





/ dr'pf(r') 


A=l 





-j*> 



2s 



ds ds 2 
dV^dpfiXr') 



ds 



dX 



- VJs] 



d 2 pf(Xr') 
dX 2 



(3.13) 



J A=l 
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Inserting these results into Eq. ( |3.10| ) for Cm we end up with 



C 



M 



dr 



~*dp s dpf 

~ m ^x +9s ^x 



dr , pf(r')s^--3(H 2 + c(l> 3 ) 



;9sP 



cS 1 dr'pf(r" 



2s^ + s 2 
ds 



d 2 V a 
ds 2 



ds 



+ 460 3 + 3cq 



dX 



A=l 



+ -j9vP J dr'p{r' 
+ ^9 P P3 J dr'p 3 (r' 



2s 



ds 
dV, 



+ s* 



ds 2 



2s 



d 2 V 

p _|_ s 2 p 



ds ds 2 

After some algebra it is possible to recast the restoring force of the monopole state as 



(3.14) 



C 



M 



dr 



+ 3 ( m s<^ + -H 



OX 



m z v V' z - m 2 p R 2 



(2m s + b(p )-— + 2m v \/— + 2m p R— 



(3.15) 



A=l 



Note that as in the case of 7^ a the same expression holds in both the RTF and RETF 
models. In each model one just has to compute Cm with the ground-state densities and fields 
obtained from the solution of the corresponding variational equations (which are modified 
by the inclusion of the corrections £2 and p^)- The derivatives with respect to A entering 
Eq. ( p.!5|) can be calculated as indicated in Appendix B, Eqs. ( |B.4j )- ([B.7| ). We have found 
Cm to take positive values for all of the (linear and non-linear) parameter sets we have used 
in the RTF and RETF calculations. A large part of the final value of Cm (usually far more 
than a half) is due to the contribution of the term —mdp B /dX\x=i- 

Evaluation of the integrand of Eq. fl3.15| ) in the limit of symmetric infinite nuclear matter 
gives the result 



Koop c 



k 2 

fir, 00 



^Pto + 3 
m 2 e FiOC 



Poo + 3- 



dm* 



-Pc 



OX 



A=l 



^F.oo 

From the field equation for the scaled scalar field in nuclear matter we have 

\ -1-1 

2 / Ps,oo 



dm* 



dX 



-3& 



A=l 



— O 

s roo 

eF,oo 



m s + 3g z s 



ml 



Poc 
^F,c 



2b<Pt a + 3c« 



00 



(3.16) 



(3.17) 



and, as expected, in ( |3. 16| ) is seen to coincide with the expression of the bulk nuclear 



incompressibility in the relativistic model [p7|,|38 
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3.3 Mass parameter 



As explained in Appendix A the mass parameter of the giant resonance is obtained from 
the second derivative of the scaled energy in a moving frame with respect to the collective 
velocity A = dX/dt, see Eq. ( |A.13| ). The relation between the collective velocity A and the 
velocity v of the moving frame is provided by the continuity equation (|A.15|) for the scaled 
system. This equation suggests a radial velocity field of the form v = —\u(r)r/r up to first 
order in A for the monopole mode [p^|13|) |15l . In terms of the displacement field u(r) the 
mass parameter ( |A.13| ) is written as 



B = J dru 2 {r)H, 



while the continuity equation (|A.15| ) becomes 

dp\(r) 



dX 



- V 



p x {r u{r - 

r. 



At A = 1 Eq. ( p,19|) determines the displacement field as 

1 



u(r) 



where 



PtO) 



p(r)r 2 Jo 

dp x (r 



dr'r' pt{t'), 



dX 



(3.18) 



(3.19) 



(3.20) 



(3.21) 



A=l 



is the so-called transition density. 

For the monopole mode p\(r) = X 3 p(Xr) and, thus, the transition density is given by 

dp(r) 



Pt(t) = 3p(r) + r- 



dr 



(3.22) 



which is known as the Tassie transition density. Partial integration of Eq. ( |3.20| ) with 



(|3.22| ) leads to the well-known result u(r) = r for the displacement field under the scaling 



transformation. The mass parameter of the monopole oscillation thus becomes 1 3 1,32 



B 



M 



drr 2 H. 



(3.23) 



Finally one calculates the excitation energy of the monopole state in the scaling model as 
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Let us mention in passing that in the non-relativistic approach the mass parameter is 
derived as B^ = J drr 2 mp, with m being the nucleon mass. Actually, provided that Bm 
is replaced by By[ in (|3.24 ), the scaling energy of the monopole vibration can be formally 



expressed as in the non-relativistic sum-rule approach. That is, = ym 3 /mi, where the 
moment mi is the energy weighted sum rule 

m a = -A(r 2 } = (3.25) 
m m 

and rriz is the plus three energy moment related to the second derivative of the scaled energy 



3.4 Constrained calculation 



The giant monopole resonance can also be studied by performing a constrained calculation. 
In the semiclassical context one has to minimize the constrained functional 



J dr[H — r]r 2 p] = E(rf) — i] J drr 2 p 



(3.26) 



with respect to arbitrary variations of the proton and neutron densities and of the meson 
fields. The densities, fields and energy obtained from the solution of the variational equations 



associated to ( 3.26|) , now depend on the value of the parameter r]. The nuclear r.m.s. radius 
is calculated as 



Rr, 



— J drr 2 p 



1/2 



(3.27) 



where A is the mass number of the nucleus. 

The constrained energy E{rj) has a minimum at rj = which corresponds to the ground- 
state r.m.s. radius Rq. Following Refs. [ 12], 13 , [14]jL5| one expands E(rj) in a harmonic ap- 
proximation about Rq to obtain the constrained incompressibility of the finite nucleus as 



K C A = ^R 2 



1 ^ 2 d 2 E{r 1 ) 
,S"° 8R 2 



(3.28) 



ri=0 



In the constrained model the displacement field is also given by Eq. ( p.20[) ; now with a 
transition density f>r(r) = dp(r,r))/ds\ v= o, where s = R v /Ro — 1 denotes the collective 
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variable of the constrained monopole oscillation |L2],[T3],[U|] . The frequency of the constrained 
isoscalar monopole vibration is computed as 



E c u = ^ • (3-29) 



Again, we may notice that if in this equation the inertia parameter is replaced by its non- 
relativistic limit B 1 ^, the energy of the constrained monopole vibration can be nominally 
written in terms of sum rules. In the present case one hase = ^/mi/m_i, with m_i 
being the inverse energy-weighted sum rule 



2 dr) 

3.5 Numerical results 



I dR 2 v 1 d 2 E( v ) 



(3.30) 

»?=o 



Our RTF and RETF results for the excitation energy of the isoscalar giant monopole res- 
onance (GMR) obtained in the scaling approach are displayed in Table 1, together with 
the empirical estimate E x ~ 80/A 1 ^ 3 MeV [j4l| . We have considered the nuclei 40 Ca, 90 Zr, 



116 Sn, 144 Sm and 208 Pb for which recent experimental data on the GMR are available [42 . 
in addition to 16 and 48 Ca. We have employed the non-linear parameter sets NL1 (incom- 



pressibility = 211 MeV, effective mass in nuclear matter m^/m = 0.57) |43[], NL3 |4[] 
= 272 MeV, rrfjm = 0.60), NL-SH || (K^ = 355 MeV, m*Jm = 0.60) and NL2 



(/Too = 399 MeV, m^/m = 0.67). The predictive power of these parametrizations is 



well known and some examples can be found, e.g., in Ref. jRJ and references quoted therein. 
The relatively new parameter set NL3 is considered to be the most successful relativistic 
effective interaction so far. It is to be noted that the RMF parameter sets are determined 
by least-squares fits to ground-state properties like radii, binding energies and spin-orbit 
splittings of a few spherical nuclei. Then, there is no further adjustment to be made in the 
parameters of the interaction. 

From Table 1 we see that the smaller the mass number, the larger is the monopole 
excitation energy. The energy of the GMR increases with increasing in the various 
parameter sets. For example, in the RETF calculation the excitation energy for 208 Pb is 
12.7 MeV with NL1, while it is 18.4 MeV with NL2. At first glance the dependence on 
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is roughly linear for each nucleus. We realize that, overall, the importance of the gradient 
corrections of order h 2 included in the RETF approach is small for the GMR energy: the 
RETF energies differ by less than 10% from the RTF energies. The largest deviations 
between the RETF and RTF results are found for the lighter systems, where the surface 
terms are comparatively more important. If we analyze the relative difference between the 
RETF and RTF energies it is seen to decrease with increasing mass number in all sets, with 
the sole exception of 16 with NL3. For example, the difference goes from —7% in 16 to 
— 1.6% in 208 Pb with NL1. We observe that the sign of the correction to the energy of the 
monopole state due to the h 2 terms depends on the value of K^. In the case of NL1 the 
change of RETF with respect to RTF is negative. For NL3 and NL-SH the change becomes 
more and more positive with (the effective mass m^/m of these two sets being almost 
the same). When we look at NL2 the change is again positive, but smaller in relative value 
than for NL-SH owing to the larger effective mass of NL2, which tends to counterbalance 
the effect of K^. 

It is usually recognized that microscopic calculations of the isoscalar GMR energy in 
nuclei provide a reliable source of information on the nuclear matter incompressibility 
47|,48| . The value of is an important ingredient not only for the description of finite 
nuclei but also for the study of heavy ion collisions, supernovae and neutron stars. In prac- 
tice one has several effective interactions which differ mostly by their prediction for Koo, but 
otherwise reproduce satisfactorily the experimental data on ground-state properties. Com- 
parison of the calculated GMR energies with experiment restricts the range of acceptable 
values for the nuclear matter incompressibility of the effective nuclear force. From Table 1 
we see that the empirical law E x ~ 80 /A 1 / 3 MeV roughly lies in between the predictions 
of the NL1 and NL3 parameter sets, as expected from the reasonable value of in these 
interactions. On the contrary, the NL-SH and NL2 parametrizations have too high a value 
of Koo and clearly overestimate the empirical curve and the experimental data for all the 
considered nuclei. In Figure 1 we have drawn further RETF results for the excitation energy 
of the monopole mode in comparison with the experimental data listed in Ref. ff9|| , as a 
function of the number of particles of the nucleus. (The RTF calculation displays basically 
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the same trend as the RETF results.) Again, it is clear that the NL-SH and NL2 sets are 
unable to describe the experimental values. The experimental points are roughly enclosed 
within the predictions of the NL1 and NL3 sets. For medium mass nuclei NL1 is closer to 
experiment than NL3, while for heavier nuclei the experimental energies tend to approach 
the results obtained with NL3. 

A further inspection of Table 1 shows that no RMF parameter set seems capable of 
reproducing the mass-number dependence of the experimental data over the whole analyzed 
region, particularly in the lighter nuclei. One should note, however, that the scaling cal- 
culation provides a prediction for the mean value or centroid of the excitation energy of 
the resonance. To establish a relation between the incompressibility and the experi- 
mentally measured energies of the monopole mode the most favourable situation is met in 
heavy nuclei, where the strength of the GMR is less fragmented than in medium and light 
nuclei |12,4~5]. If we take into account the excitation energies of 116 Sn, 144 Sm and 208 Pb, 
according to Table 1 the nuclear matter incompressibility of the relativistic interaction 
should lie in the range 220-260 MeV. If we only consider 144 Sm and 208 Pb, as in Ref. f20|j , 
then would be restricted to the range 230-260 MeV. For comparison, the authors of 
Refs. |TT1 , |20[] conclude that the value of should be close to 250-270 MeV from their 
time-dependent RMF and relativistic RPA calculations. The analysis of the relativistic 
RPA peak energies reported in Ref. || for 208 Pb suggests instead a range 235-250 MeV for 
Kqq. On the other hand, non- relativistic Hartree-Fock plus RPA analyses using Skyrme 
and Gogny interactions determine to be 215 ± 15 MeV [47,4^], thus lower than in the 
RMF model. 

The restoring force Cm of the GMR defines the incompressibility K S A of the finite nucleus 
in the scaling model through Cm = AK S A . In the limit of an arbitrarily large nucleus 
K A should approach the nuclear matter value K^, see Eq. ( |3.16| ). This suggests a linear 
dependence of the incompressibility of finite nuclei on the bulk incompressibility of 
the effective interaction. In Figure 2 we display the value of K A for the nuclei 16 0, 40 Ca, 
116 Sn and 208 Pb obtained from our RETF calculation, as a function of K^. Apart from the 
parameter sets discussed in Table 1, we have employed the sets NL-Z, LZ, LI, L2, L3 and 
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HS compiled in Table 3 of Ref. [60L the sets NLB, NLC and NLD from Ref. [51 , NL-RA1 



from Ref. |52j and LO from Ref. |43]. The bulk incompressibility of these sets spans a wide 



range of values, from ~ 175 to 625 MeV. The sets LO, LZ, LI, L2, L3 and HS correspond 
to the linear model without scalar self-interactions. 

The results of Figure 2 show a remarkable linear behaviour of K S A with the compression 
modulus Krx,. This is more true for the heavier systems on the one hand, and for the non- 
linear parameter sets on the other hand. The linear sets show a considerable dispersion, 
but one should take into account that only LO and LZ were optimally adjusted to nuclear 
ground-state properties and that, furthermore, LI, L2 and L3 do not include the p field. 
The incompressibilities of 208 Pb and 116 Sn are nearly the same. To see a perceptible change 
one has to go to 40 Ca. In Figure 3 we have drawn the excitation energy E^ of the monopole 
state versus the incompressibility. As expected from the pattern displayed by K\, the 
monopole excitation energy increases smoothly with the bulk compression modulus, roughly 
as a linear function of the square root of (in agreement, e.g., with Refs. PJ47|). Both 
the effective mass at saturation m*^ and the mass of the scalar meson m s play a major role 
in the determination of the nuclear structure properties in the RMF theory. The effective 
mass has a direct influence on the spin-orbit force and the single-particle levels, while the 
scalar mass is related with the range of the attractive part of the effective nuclear force and 
thus strongly affects the nuclear surface. One may wonder whether these two quantities 
have some effect on the energy of the breathing mode. Figure 4 shows that this is not the 
case, as no evident correlation seems to exist between and the value of or m 8 . 

From the data represented in Figure 2 we obtain the relations 



n A 


= 0.66^ 


- 12 MeV 


for 208 Pb, 






= 0.65^ 


-9 MeV 


for 116 Sn, 




n A 


= 0.57^ 


-7 MeV 


for 40 Ca, 






= 0.45^ 


-8 MeV 


for 16 0. 


(3.31) 


gen 


and calcium 


are not as 


meaningful as 


for the heavier nuclei, 



and we give them mostly for the purpose of illustration. Though there is a dependence on 



the mass number, the slope of the linear fits ( 3.31 ) is visibly smaller than unity. Moreover, 
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we have obtained a non-vanishing constant term. This is consistent with the leptodermous 
expansion of the incompressibility of a finite nucleus, inspired from the liquid drop for- 
mula, in which one separates the volume, surface, symmetry, Coulomb and higher-order 
contributions by writing 

K\ = K 00 + K s ami /AV 3 + K sym (N - Zf/A 2 + K Coul Z 2 /A 4 / 3 + ■■■. (3.32) 

The total incompressibility K\ receives a sizeable contribution from the surface term and 
smaller contributions from the symmetry and Coulomb terms The sign of these 

terms is negative and they considerably decrease the value of K S A in real nuclei with respect 



to the i^oo limit ||17| , |18| , |53| , |54f . A key point is the fact that the surface incompressibility K^ ur{ 
in the scaling model varies almost as a linear function of K^, which guarantees that the 
surface effects do not destroy the regular behaviour of K\ with K^. For instance, our RETF 
calculations of the K* ur{ coefficient for several relativistic parameter sets [pH[j , by using the 
scaling method in semi-infinite nuclear matter, confirm that i^| ur f in the relativistic model 
indeed behaves roughly linearly with the bulk compression modulus, as happens with non- 
relativistic Skyrme forces. In fact we have found |55| that the rule of thumb K^ nrf ~ — 



known from the non- relativistic approach ||54|| , also applies to non-linear RMF parameter 
sets having not too large values of the compression modulus. In the case of the linear a — uj 
sets, which have a high value, one instead finds K^ mf ~ — l^K^. 

It is interesting to compare (|3.31|) with the equations K S A = 0.627^00 + 23 MeV for 208 Pb 



and K A = 0.49/0,0 + 35 MeV for 40 Ca obtained in Ref. |56| from an analysis of the scaling 



incompressibility performed with the Skyrme forces SkM*, SGI and SIII. In the relativistic 
model the independent term of the linear fits is negative and seems to be more constant 
with mass number, but the coefficient in front of Koc is similar in both the relativistic and 
non-relativistic model. The authors of Ref. signaled that the slope obtained for 208 Pb 



with the Skyrme forces approaches the hydrodynamical value 7r 2 /15 = 0.658, though they 
stressed that this might be just accidental. It is at least curious to come across with the 
same value in the relativistic model. 

In assuming a nuclear matter approach Nishizaki et al. estimated the monopole 
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excitation energy of a finite nucleus in the relativistic model as 



^ = V^y (3 - 33) 

where /Xqq is the chemical potential of nuclear matter (including the nucleon rest mass), 
(r 2 ) = |i? 2 and R = 1.2A 1 ' 3 fm. They evaluated ( |3.33D for the linear model of Walecka 
using /loo = 923 MeV and = 525 MeV and found Efa = 160/A 1 / 3 MeV. This result has 
the correct dependence on the mass number but it is twice as large as the empirical value 
~ 80 /A 1 / 3 MeV. Calculating ( ggp for the non-linear parametrizations NL1, NL3, NL-SH 



and NL2 one finds = 102, 115, 132, and 140/A 1 / 3 MeV, respectively. Compared to the 
empirical law these values are too large by a factor ~ 1.5-1.8, depending on the compression 
modulus of the force. If we furthermore compare with the results of the calculations for 
finite nuclei listed in Table 1, we realize that the finite size effects reduce the prediction 
obtained from nuclear matter by a factor ranging from 1.6 in 16 to 1.3 in 208 Pb, almost 
independently of the parameter set. 

In a recent work Piekarewicz |§ has given a thorough presentation of the relativistic RPA 
formalism and has computed the isoscalar monopole mode for several closed-shell nuclei. In 
the numerical calculations he has used the non-linear sets NLC and NLB and the linear 
set L2', which have the nuclear matter incompressibilities = 224, 421 and 547 MeV, 
respectively. We present in Table 2 our values obtained from the scaling method versus the 
RRPA peak energies of the isoscalar mode taken from Ref. ||, for the systems 40 Ca, 90 Zr 
and 208 Pb. A fairly good agreement is found between our semiclassical calculations and 
the more fundamental RRPA approach. The differences are well below 5% in 208 Pb and, 
excluding the RETF result for 40 Ca with L2', below 10% in 40 Ca and 90 Zr. As discussed in 
Ref. 0, it becomes difficult to even identify a genuine resonance in the RRPA distribution 
of the isoscalar monopole strength for medium-size nuclei such as 40 Ca with the parameter 
sets NLB and L2' which have large compression moduli. 

The GMR has also been studied by means of constrained calculations in the RMF model 
and, based upon them, with the more ellaborate generator coordinate method (GCM). The 
constrained calculations in our semiclassical approach (see Section |3.4| ) are carried out in a 
similar way to that of Refs. [12 ,1^,14 ,U| within the quantal Hartree approach. We report in 
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Table 3 the excitation energies of 16 0, 40 Ca, 90 Zr and 208 Pb calculated with the constrained 
RETF method (for the NL1, NL2, NL3 and NL-SH sets), besides the constrained RMF 
(CRMF) Hartree results of Ref. [H] and the GCM results of Ref. [ Il|. In non-relativistic 
RPA calculations it is common to utilize the moments of the strength function to analyze 
the monopole resonance |L9||. The lowest moments correspond to simple sum rules and in 



the limit of small amplitude oscillations the ratios Jm^/mi and Jmi/m-i can be identified 



respectively, with the scaling and constrained monopole excitation energies fli8| , [47| . |57 . 

We see that the excitation energies of the monopole state are smaller in the constrained 
model (Table 3) than in the scaling model (Table 1). This is in agreement with the non- 



relativistic RPA inequality ymi/m_i < ym 3 /mi [[HJ. When the comparison is possible, 
the energies obtained with the GCM are systematically smaller than in the CRMF Hartree 
model, which in turn are smaller than in the constrained RETF approach. The constrained 
RETF values agree very well with the GCM and CRMF values for NL1, but the agreement 
worsens for light nuclei with the other parameter sets. In the case of 208 Pb, for which the 
semiclassical technique should work better, the constrained RETF calculation overestimates 
the GCM value by around 1 MeV, the same magnitude by which the CRMF and GCM results 
differ (for NL1 and NL-SH). 

Vretenar et al. Jll| have studied the GMR with the time-dependent RMF approach. (We 
recall that very recently it has been demonstrated that the relativistic RPA, with inclusion 
of Dirac sea states, is equivalent to the small amplitude limit of the time-dependent RMF 



theory in the no-sea approximation 0].) In the calculations of Ref. [J 1 the main peak 
appearing in the monopole strength distribution of 208 Pb has energies 12.4 (NL1), 14.1 
(NL3), 16.1 (NL-SH) and 17.8 MeV (NL2), while in the case of 90 Zr the energies are 15.7 
(NL1) and ~ 18 MeV (NL3). For 208 Pb our scaling results (cf. Table 1) show a good 
agreement in all the parameter sets. In fact, if we focus on the RETF values, we see that 
the scaling energies are an upper bound of the time-dependent RMF energies, while the 
constrained energies of Table 3 represent a lower bound (apart from the case of NL-SH by a 
little deviation). For 90 Zr, however, both the scaling and constrained semiclassical excitation 



energies are larger than those of Ref. ||11|| . It should be pointed out that the Fourier spectrum 
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of Zr in the time-dependent RMF calculation is considerably fragmented (specially for the 
sets with higher K^) and then the determination of the centroid energy remains more 



uncertain 11 
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4 Giant quadrupole resonance 



In the quadrupole vibration the particle density scales as [19 



p x (r) = p(x/\,y/\,\ 2 z). (4.1) 

While the volume element is conserved in both coordinate and momentum space, the mo- 
mentum distribution, which remained spherically symmetric in the monopole oscillation, 
becomes highly deformed in the quadrupole case PJSBf: 

Px = {Xp x ,Xp y , Pz /X 2 ). (4.2) 

One has to note that the spherically averaged form of the distribution function TZ(r, p) 
cannot be employed in the quadrupole scaling calculations due to the deformation of the 
Fermi sphere [j58|| . This means, in particular, that the spherically symmetric expressions 
(|2.4| ), (|2.5|) and ( [2.12|) of the semiclassical energy density and scalar density are no longer 
valid for use in the quadrupole scaling. That is, first, one should replace p by p x in the 
semiclassical expansion of the relativistic distribution function TZ(r,p) |M and then obtain 



its moments (energy and densities) as a function of the collective coordinate A. This ex- 
traordinarily complicates the expressions if the distribution function with terms up to order 
h 2 is to be used. Since the final magnitude of the contribution of the 7i 2 -order corrections in 
the semiclassical calculation of the excitation energy of giant resonances is not very signifi- 
cant, we will work at the Thomas-Fermi level in the present study of the giant quadrupole 
resonance. 

In the Thomas-Fermi approach the relativistic distribution function is proportional to a 
step function (Appendix A and Ref. which vanishes for single-particle energies above 



the Fermi level. The Thomas-Fermi energy density of the non-linear a — u model after 
scaling then reads 



Kx = E / dP\/pl + m *x 2 © (/V - \fp\ + ml 



it. 



+ \g^\P e i + \b<Pl + ^c0l + ~£ V V A PA + \g P R\pi\ + \eA x p pX , (4.3) 
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where denotes the step function, fj, q \ is the chemical potential of the scaled system for 
each kind of nucleon and the single-particle potential u q \ is given by 



u q \ = g v V x + -g P R\T 3 + -eA\(l + r 3 ). 
The scaled effective scalar density plf has been defined through 
9sPsX - b4>\ - C(p\ 



(4.4) 



9spfx 



9sm x 



6 (u qX - y/ft + mf - Uqx^j - b<p 2 x - c<p\. (4.5) 



(2vr 7 „ 

The position and momentum variables in these expressions scale according to the rules ( [4.1| 



and ( fj~2[) in the quadrupole case. 

To obtain the restoring force Cq of the quadrupole oscillation we have to compute the 
second derivative of the scaled energy with respect to the collective coordinate A. The first 
derivative reads 




8J 



(2*) 



dr 

+ Tx ldr 



dpx 



Vx + m *x 2 ^ 



© [PqX 



PX + m f - UqX 



9sP s x dx 



1 



1 



1 



+ 7;9vV x px + ~9 P R\P3\ + ~eA\Pp\ 



(4.6) 



It can be checked that this equation identically vanishes at A = 1, as in the non-relativistic 

case 



19|| . Before deriving again ( [4.6D it is helpful to take into account that, for instance, 



drpx(r)Vx(r) = J drp(^,A 2 z) J dr 'p(j^ 
drp(r) J dr' P (r')VUs\), 



X z' 



4/T \r 



where 



VMx) 



9v e 



-m v s x 



sx = (Ax — Ax', \y — Xy', zjX — z'/X ). 



(4.7) 



(4i 



4tt s a 

With this, after some algebra, the restoring force of the quadrupole mode can be put in the 
form 

d 2 r 1 
- / drHx(r) 

J A=l 
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s ds 



(4.9) 



After performing the angular average in the integral over p and in the integrals over r 
and r', we finally get 
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(4.10) 



As far as the nuclear part is concerned this result coincides with the one derived in Ref. [[31 



for nuclear matter using a local Lorentz boost and the scaling method. The contributions 
from the meson fields agree with the result obtained from the potential part of an effective 
density-independent nuclear force in the non-relativistic model ]57j (and the contribution 



from the Coulomb field agrees with that given in Ref. [ 19] ) . 

As in the monopole oscillation to calculate the mass parameter one needs the continuity 
equation in a moving frame, Eq. ([A.15Q . For the quadrupole vibration we have p\(r) = 
p(x/X, y/X, \ 2 z) and the continuity equation (|A.15|) is fulfilled by v = —X(—x,—y,2z) = 



— AV[y 47r/5r 2 Y"2o(^)] at A = 1 [31], which provides the connection between the velocity 
of the moving frame and the collective coordinate. Proceeding similarly to the monopole 
case, i.e., inserting this velocity field into Eq. (|A.13 ) and taking the second derivative with 
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respect to A, the mass parameter of the quadrupole mode is found to be 

B Q = 2 / drr 2 H, (4.11) 



assuming the nucleus to be spherical. The excitation energy of the quadrupole state then is 



£q = ^- ( 4 - 12 ) 



The transition density in the quadrupole case is given by 
dp\{r 



= V PA (r)| A=1 ■ V[ v ^75r 2 F 20 (fi)] = J^r^Y^n), (4.13) 

A— 1 



where again we have assumed the density to be spherically symmetric at A = 1. 
4.1 Numerical results 

As we have indicated, our calculations for the quadrupole mode are restricted to the RTF 
approximation. We collect in Table 4 the calculated excitation energy of the quadrupole 
oscillation for 16 0, 40 Ca, 48 Ca, 90 Zr and 208 Pb, along with the empirical law E x ~ 



MeV and some experimental data taken from Ref. J4TJ. The theoretical results shown in 
this table correspond to the non-linear sets NL1, NL3, NL-SH and NL2, and to the set LZ 
(Koo = 586 MeV, mj^/m = 0.53) which we take as a representative of the linear sets. 

One can see that the four non-linear o — uj parametrizations reproduce the empirical 
trend and that, contrary to the situation found in the monopole case, they give rather 
similar results for each nucleus. This is due to the fact that the energy of the quadrupole 
vibration is basically independent of the bulk compression modulus of the effective force. 
Nevertheless, the comparison with experiment favours the NL3 set and, especially, the NL1 
set (i.e., those sets with a lower incompressibility) . In fact, if the incompressiblity of the 
force is very large (set LZ) the theoretical predictions clearly overestimate the experimental 
values. The relativistic results of the non-linear sets compare well with those obtained 
in non-relativistic Hartree-Fock and extended Thomas-Fermi calculations using Skyrme 



forces [26]. Calculations of the isoscalar giant quadrupole resonance are rather scarce in the 



relativistic domain. Time-dependent RMF calculations of this mode have been carried out 
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in Ref. []X0[| using the NL-SH parameter set. Our relativistic Thomas-Fermi calculation is in 
good agreement with the excitation energies of 23.6, 17.7 and 17.7 MeV for 16 0, 40 Ca and 
48 Ca, respectively, reported in that work. 

The energy of the quadrupole excitation has also been evaluated by Nishizaki et al. [[H 
from a nuclear matter approach as 



E 8 r 



6*Lo 



Q = \ i ^r, (4.14) 

where (r 2 ) has been defined in Eq. ( p.33| ). In this approximation the restoring force of the 
quadrupole vibration corresponds to the nuclear matter limit of Eq. (|4.10 ), where all the 
terms with derivatives of the meson fields vanish and only the first term survives. Note that 
the incompressibility of the interaction does not enter Eq. (|4.14 ). According to this 
equation one obtains E S Q = 85, 84, 81, 80, and 76 /A 1 / 3 MeV for the LZ, NL1, NL3, NL-SH 
and NL2 sets, respectively. We thus see that in nuclear matter .Eq decreases as the value 
of the effective mass at saturation of the force (m^) increases. However, in the full RTF 
calculation for finite nuclei (Table 4) the regular pattern of Eq with observed in nuclear 
matter is destroyed by the finite size effects. In the case of finite systems one not only has 
the additional contribution from the meson fields into Eq. ( 4.10|) , but also the nuclear part 
is modified by the shape of the nuclear surface, this one depending in turn on the mass of 
the sigma meson m s . Such effects mask the simple relation of Eq with shown by the 
naive nuclear matter approximation. 
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5 Summary and conclusions 

We have studied the isoscalar giant monopole and quadrupole resonances of finite nuclei 
by means of the scaling method and the Thomas-Fermi and extended Thomas-Fermi ap- 
proaches to relativistic mean field theory. Self-consistent numerical calculations for realistic 
non-linear a — uj parameter sets have been discussed. Previous relativistic investigations 
with the scaling method either relied on a leptodermous expansion of the finite nucleus 
incompressibility [ p!3|JTB|JT7[ ] , or were limited to the linear o — uo model for symmetric and 



uncharged nuclei at the Thomas-Fermi level [31,32 



In the present approach one starts by scaling the spatial and momentum coordinates 
of the semiclassical distribution function in a moving frame. By taking the derivatives of 
the scaled energy in the moving frame with respect to the collective coordinate and the 
collective velocity, one obtains the expressions from which the restoring force and the mass 
parameter of the resonance can be computed. The underlying reason for the success of the 
method is that in the semiclassical approach the energy functional is written explicitly in 
terms of the local Fermi momentum and of the local effective mass, which allows one to 
easily perform the scaling. Due to the finite range of the relativistic interaction no compact 
formulas can be obtained as in the case of non-relativistic Skyrme forces. Nevertheless, the 
scaling excitation energies of the monopole and quadrupole resonances only depend on the 
ground-state densities and fields, which means that they can be computed as a by-product 
of a semiclassical self-consistent calculation of the ground state. 

We have found that the total contribution to the excitation energy of the GMR coming 
from the gradient corrections of order h 2 , which are included in the RETF approach, does 
not modify the Thomas-Fermi result very much. The strength and sign of these corrections 
of order h 2 is strongly correlated with the nuclear matter incompressibility and the effective 
mass at saturation of the relativistic interaction. 

We have investigated the relation between the incompressibility K S A of finite nuclei in 
the scaling model and the compression modulus of nuclear matter K^, employing a vari- 
ety of relativistic parameter sets. The dependence is roughly linear, as in non-relativistic 
analyses. Even a nucleus such as 208 Pb is not large enough to obtain a relation of propor- 
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tionality between K\ and K^. The excitation energy of the monopole oscillation increases 
smoothly with K^ 2 , in correspondence with the behaviour of K\. No regular pattern of the 
monopole excitation energy with the mass of the scalar meson or with the effective mass of 
the interaction has been observed. 

The experimental excitation energies of the monopole oscillation in medium and heavy 
nuclei lie in between the results obtained with the NL1 and NL3 parameter sets. An analysis 
of the calculated breathing-mode energies for 116 Sn, 144 Sm and 208 Pb, for which precise 
experimental data exist, predicts that the nuclear matter incompressibility should be around 
220-260 MeV (230-260 MeV if only 144 Sm and 208 Pb are taken into account). A similar 
analysis carried out in Refs. |]TT| , |2"0"|] using time-dependent RMF and relativistic RPA results 
predicts a value slightly higher: 250-270 MeV. From the relativistic RPA peak energies 
given in Ref. for 208 Pb we extract a range of 235-250 MeV. Thus, all these relativistic 
calculations point to a value of roughly 250 ±20 MeV for K^, which is higher than the non- 



relativistic estimate of 215 ± 15 MeV from Skyrme and Gogny forces [pP7| , ^q| . Relativistic 
parameter sets with large values of (such as NL-SH or NL2), which may otherwise 
perform well in describing the data for nuclear masses and radii, should be discarded on the 
basis of the experimental information on breathing-mode energies. 

The results computed with the scaling method represent an upper bound of the mean 
excitation energy of the GMR, to the extent that they are related with the cubic weighted 
sum rule. Instead, the breathing-mode energies obtained from constrained calculations 
rather represent a lower bound, since they are related with the inverse energy-weighted sum 
rule. Actually, with all the parameter sets and nuclei analyzed, we have found the calculated 
monopole energies to be larger in the scaling approach than in the constrained approach. 

Our calculations of the excitation energy of the quadrupole oscillation have been re- 
stricted to the Thomas-Fermi approach, to simplify the problems related with the distor- 
tion of the Fermi sphere. All the considered non-linear parameter sets reproduce fairly well 
the empirical trend, rather independently of the value of the compression modulus of the 
force. Although a nuclear matter estimate predicts a decrease of the quadrupole excitation 
energy with an increase in the value of the effective mass at saturation, the finite size effects 
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and additional contributions from the meson fields mask this trend in the self-consistent 
calculations for actual nuclei. 

In conclusion, we hope to have shown that the scaling method can be confidently used to- 
gether with the relativistic Thomas-Fermi approach to estimate the excitation energy of the 
isoscalar monopole and quadrupole resonances in a simple and reliable way. The results for 
the breathing mode turn out to be in good agreement with the outcome of dynamical time- 
dependent RMF and relativistic RPA calculations. We can thus conclude that, similarly 
to the non-relativistic case, also in the relativistic framework the semiclassical excitation 
energies obtained with the scaling method simulate the results of the RPA. The method 
introduced in this work also allows one to self-consistently compute the surface incompress- 



ibility coefficient for relativistic interactions [pq| . The study of other multipolarities using 
a generalized scaling simultaneously with the relativistic Thomas-Fermi approach may be 
a worthwhile task to pursue. 
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Appendix A 



In this appendix we derive the Thomas-Fermi expression of the energy of a nucleus described 
by the non-linear a — u model in a frame moving with velocity —v. As a final product 
we obtain general equations for the restoring force and the mass parameter of the giant 
resonance. For simplicity we shall consider an uncharged symmetric nucleus (the p meson 
field and the electromagnetic field behave like the vector field), and shall not include the 
corrections of order h 2 to the Thomas-Fermi approximation. 

The semiclassical expressions of densities and energies are most conveniently derived from 
the so-called phase-space distribution function ||. For a Hamiltonian a ■ p + (3m* + g v V 
the distribution function in Thomas-Fermi approximation reads [3~3] 



n 



le(v-e-g v V) 



I + -{ct ■ p + f3m*} 
e 



(A.l) 



where \i is the chemical potential, e = yp 2 + rn* 2 and I is the 4x4 unit matrix. Due to the 
step function in ( |A.1| ), p takes values from zero up to the Fermi momentum pp. The scalar 
field {(f)) and the time-like component of the vector field (V) transform to a frame which 
moves with velocity — v like <fi' = <fi and V = 'yV, with 7 = 1/y/l — v 2 . The distribution 
function in the moving frame then is given by 



K' = \Q{p' - e' - 19v V) [/ + ^{a • (p' - g Y V) + [3m*} 



(A.2) 



where // is the chemical potential in the new frame, e' 
defined 

V' = 7 vV. 

It is easy to see that 9(/i' — e' — 75^) = 0(/x — 
and fi' through their values in the rest frame: 



p' — g v V') 2 + m* 2 , and we have 

(A.3) 

6(j»F — p)] by expressing e' 



e' = 7(e +p ■ v), // = 7(// + p ■ v). 



(A.4) 



The baryon density in the moving frame is obtained as 
dp' 



P 
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(A.5) 



where we have taken into account that dp' /e' = dp/e (Lorentz scalars) and the fact that the 
trace of the distribution function 1Z' equals 20(pp ~ p)- Similarly, the transformed scalar 
density is 



Ps 



dp' 



Tr[pn'] 



dp m* 



(2tt) 3 l ' " j J (2it) a e 
The energy density in the moving frame is given by 

dp' 



(A.6) 



w 



(27T) 



-Ti\H'n' 



V'ch) 2 + mid 2 



7 2 ( Wf + ^mlV 2 - (V x V') 2 - m 2 v V' 2 

where H' = a ■ (p' — g v V) + (3m* + jg v V. 

If spherical symmetry of the meson fields is assumed Eq. QA.7 ) becomes 

dp 7 2 



(A.7) 



(2vr) 3 e 



[e + p ■ v)(e + p ■ v + g v V)&(p F - p) 



1 

+ 2 

1 

~ 2 



:V0) 2 + m s 2 2 + 37V(V0) 2 



+ |&^ 3 + i C 4 



:VF) 2 + m 2 \Z 2 + - 7 V(V\/) 2 
3 



(A.8) 



After integration over momentum, the relativistic energy density in the moving frame can 
be written as 



n' = 7 2 In + v 



pe F + g vP V + ^(V0) 2 - i( VV) 2 - W 



(A.9) 



where ?^ is the energy density in the rest frame (Section 0). Equation ( |A.9 ) agrees with 
the transformation law of the stress tensor as discussed in Ref. |59| . For a uniform system 
(V0 = VV = 0) it also coincides with the result obtained in Ref. |3l| from a local Lorentz 
boost. 



Finally, the energy of the system in the moving frame is obtained by integrating 
over the space. Taking into account the Lorentz contraction of the volume element, this 
yields 



E(v) = J ^-n' = J d ri |(1 - v 2 )H + v 2 



pe F + s vP V + ^(V0) 2 -^(VVf 



. (A.10) 



31 



Combining this result with the meson field equations and the virial theorem derived in 
Appendix B, Eq. ( B.10|) , the energy in the new frame reads 



E(v) = J dr 7 {(l - v 2 )H + v 2 H} = J dr^H. 



(A.ll) 



The restoring force of the monopole and quadrupole oscillations is obtained by appro- 
priately scaling the densities and mean fields in Eq. ( |A.11| ) and then computing the second 
derivative at v — and A = 1: 



C 



dX 2 



f-K 

J 7 



dX 2 



drH> 



(A.12) 

j 1 ,=0,A=l L^' y " J A=l 

where 7i' x and 7i\ denote the scaled energy densities in the moving and rest frames, re- 
spectively. The mass or inertia parameter of the giant resonance is furnished by the second 
derivative of the scaled energy in the moving frame with respect to A = dX/dt: 



B 



d 2 r dr 



dX 2 



1 



OX 2 



dvyli, 



(A.13) 



J A=0,A=1 



A=0,A=1 

To evaluate ( |A.13| ) it is necessary to relate the velocity v of the moving frame with the 
collective velocity A. This is achieved by scaling the continuity equation 







dp' 



dtJ 7(2tt) 3 

which after some algebra results into 

dp\ 
dt 



V • 



dp' 



7(2?r) 



-Trfa^'l = 0, 



(A.14) 



V • (v Px ) = 0, 



(A.15) 



in terms of the scaled baryon density p\. Once the scaling law of the baryon density with 
the A parameter is specified, Eq. (|A.15| ) will provide the connection between the velocity v 
and A. 
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Appendix B 



The virial theorem results from homogeneity properties of the kinetic energy and potential 
energy components of the energy with respect to a scaling transformation that preserves 
the normalization ||60|| . For example, the scaling method has been employed to derive the 



virial theorem for the Skyrme interaction ||19| , or for relativistic particles bound in vector 
and scalar potentials [|6T|. Concerning the relativistic model discussed in the present work, 
we have given the expression of the first derivative of the scaled energy with respect to the 
scaling parameter A in Eq. (|3.8| ) of Section |3~2| . It must vanish at A = 1 (virial theorem): 

Hx(r) 







d_ 
dX 



d(Xr) 



A 3 



(B.l) 



A=l 



To evaluate the above equation knowledge of the derivatives of the scaled fields with 
respect to A is required. Starting with the omega field V\, it fulfils the Klein-Gordon 
equation 

(A-ml)V x (r) = -g vPx (r), (B.2) 



whose solution is 



Vx(r) 



e -m v \r-r'\ 



d(Xr')p{Xr')V u (s) 



(B.3) 



We employ the notation V w (s) = g v exp (— m v s)/47rs, with s = \r — r'\, as in the main text. 
On defining u = Xr and u' = Xr' one obtains V\(r) = J du 1 piu^V^u — u'|/A), whence 



dVx(r) 



dX 



A=l 



dr'p(r') s 



ds 



(B.4) 



in agreement with the result given in Ref. [57]. Analogous results are found for the scaled 
rho and Coulomb fields. The result for the scalar field is more complicated because an 
additional term appears due to the fact that the density p| ff itself is a function of A: 



dX 



Jdr'pf{r')s d ^ + Jdr'V a (s) 



A=l 



dpf(Xr') 
dX 



(B.5) 



A=l 



Since g s pf = g s p s - b4>l/X 3 - c^l/X 3 , we have 



dX 



A=l 



dps 

OX 



+ 3(60 2 + c0 3 ) - (260 + 3c0 : 



2n d<P\ 



A=l 



dX 



(B.6) 



A=l 
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with 



dps 
dX 



A=l 



5p s dm* 
5m* dX 



A=l 



5m* 



m + g s 



dX 



(B.7) 



A=l 



cf. Eq. (|3J) for dm*/dX. 

From substitution into Eq. (p.l| ) of the derivatives (|B.4|) and ( P-5| ) and of the corre- 
sponding results for the rho and Coulomb fields, on account of Eq. (|3.8|), one obtains 



= J dr 8- m*p s + -g s pf J dr 



'pf[r)s 



'\ dV<r , ,3 3 4 



1 f i i / a dVuj 1 fit / /\ dV p 1 , 
- TjQvP J dr p{r )s— - -g p p 3 J dr p 3 {r )s— + -eAp^ 



(B.8) 



Now, using for example the relation sdV^/ds = — — m Y sV W) it can be verified that 



-g vP V + m\V 2 



(B.9) 



After similar straightforward manipulations with the other fields, the virial theorem for the 
non-linear a — u model finally becomes 







dr 



£ - m*p s - -g s (j)p s - ml4> 2 - hxp* - ^c0 4 



+ \g v Vp + mlV 2 + -g p Rp 3 + m 2 p R 2 + Ka Pp 



;b.io) 



One may notice that the quantity £ — m*p s corresponds to the semiclassical average of 
Y^i l ~p\ct ■ V<y9j. Actually, in terms of the kinetic energy density r (namely, the semiclassical 
counterpart of J2i fl [ a ' V + (3m — m]ipi) we can write £ — m*p s = r + mp — mp s , which 
makes more obvious the kinetic energy component in the virial theorem. In the limit of 
symmetric infinite nuclear matter Eq. ( |B.10| ) goes over 



c * 

£•0,00 ^oo^s,oo 



2 Q s 

4 ^v„2 _ 3p = 



2mr s '°° 



- • .3 3 



(B.11) 



with P being the pressure, if equilibrium quantities are used. 

Taking advantage of Eq. ( |B.10| ) to eliminate / dr£ from the expression of J dr?i, the 
energy of a nucleus in the RMF model takes the remarkably simple form 



J dr[H — mp] = J dr 



(B.12) 
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where we have subtracted the nucleoli rest mass contribution. This expression displays 
very clearly the relativistic mechanism for nuclear binding. It stems from the cancellation 
between the scalar and vector potentials and from the difference between the scalar density 
and the baryon density (i.e., from the small components of the wave functions). We have 
verified that Eqs. ( p.lOj) and (|B.12| ) are satisfied not only by the Thomas-Fermi solutions, 
but also by the ground-state densities and meson fields obtained from a quantal Hartree 
calculation. Of course, the energy stationarity condition of the RMF model against dilation 
must be fulfilled by any approximation scheme utilized to solve the problem. 
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Figure captions 



Figure 1. The breathing-mode energies from the RETF scaling calculations are compared 
with the empirical law 80 /A 1 / 3 MeV and the experimental data reported in Ref. f49fl , 
ClS db function of the size of the nucleus. 

Figure 2. Scaling incompressibility of some finite nuclei as obtained in the RETF calcu- 
lations versus the nuclear matter incompressibility for various relativistic parameter 
sets. The value of of each set is listed in MeV. The straight lines are linear fits. 
The fit for 208 Pb is drawn by a solid line. 

Figure 3. Monopole excitation energy from RETF scaling calculations versus the nuclear 
matter incompressibility for various relativistic parameter sets. The value of of 
each set is given in Figure 2. The dashed lines are linear fits to the square root of K^. 

Figure 4. Monopole excitation energy of 208 Pb from RETF scaling calculations, as a func- 
tion of the mass of the scalar meson and of the nuclear matter effective mass of the 
relativistic interaction. The dashed lines show the sense of increasing Koo. 
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Table 1: Excitation energy of the monopole state (in MeV) obtained in the scaling approach 
by using various relativistic parameter sets (in order of increasing value of the compression 
modulus ifoo)- The experimental centroid energies are from Ref. |42[| . 





NL1 


NL3 


NL-SH 


NL2 








RTF 


RETF 


RTF 


RETF 


RTF 


RETF 


RTF 


RETF 




Exp. 


16Q 


25.1 


23.3 
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27.8 


30.7 


33.1 


34.4 


35.6 


31.7 




40 Ca 
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20.6 
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30.3 


23.4 


19.2 ±0.4 


48 Ca 
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90 Zr 
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17.9 ±0.2 
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15.9 


15.6 
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16.1 ±0.1 
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14.9 
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16.6 
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21.1 


15.3 


15.4 ±0.3 


208 pu 


12.9 


12.7 


14.5 


14.6 


16.6 


17.0 


18.1 


18.4 


13.5 


14.2 ±0.3 



Table 2: Comparison of the giant monopole resonance energies (in MeV) obtained in the 
scaling model with those obtained in the relativistic RPA ||. 







NLC 






NLB 






L2' 






RTF 


RETF 


RRPA 


RTF 


RETF 


RRPA 


RTF 


RETF 


RRPA 


40 Ca 


22.5 


22.4 


21.0 


27.7 


29.4 


27.9 


29.1 


33.0 


27.3 


90 Zr 


18.1 


18.1 


16.9 


23.3 


24.2 


24.1 


25.2 


27.4 


26.5 


208 Pb 


13.6 


13.5 


13.1 


18.0 


18.5 


18.1 


19.9 


21.0 


20.1 



40 



Table 3: Monopole excitation energy (in MeV) obtained by constrained calculations with 
various parameter sets. The constrained RMF results are from Ref. [14j] and the generator 
coordinate method results are from Ref. ITT . 



NLl NL3 NL-SH NL2 

RETF CRMF GCM RETF GCM RETF CRMF GCM RETF GCM 

6 21.8 20.9 20.2 26.0 22.6 30.0 25.8 25.0 32.4 27.1 

Ca 19.8 19.2 16.6 23.2 19.6 26.9 23.9 22.0 29.0 24.4 

3 Zr 16.5 16.3 14.1 19.1 16.9 22.1 21.1 19.5 23.7 21.9 

Pb 12.1 12.2 11.0 14.0 13.0 16.2 16.1 15.0 17.4 16.6 



Table 4: Excitation energy of the quadrupole vibration (in MeV) obtained in the scaling 



approach. The experimental values are from Ref. |41 
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11.0 


10.5 
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